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We studied the ground and excited states of 12 C based on variational calculations after spin- 
parity projection in a framework of antisymmetrized molecular dynamics(AMD). The calculations 
systematically reproduce various experimental data. It was found that the sub-shell closure and 
S'L r (3)-hmit 3a cluster components are contained in the ground state, while various 3a cluster 
structures develop in the excited states. We discussed effects of a breaking and show the importance 
of coexistence of the cluster and shell-model-like aspects. 



I. INTRODUCTION 



o 
o 

' The formation of mean-field is an important aspect in nuclear systems. It results in the shell-model-like nature of 
nuclear structure. On the other hand, in a region of light nuclei, it is well known that clustering is another essential 
feature. These two different kinds of nature, cluster and shell-model-like aspects, often coexist, compete and interplay 
between each other. The coexistence of two aspects provides a variety of structure in nuclear systems, and can be a 
key to reveal properties of light nuclei. 

12 C is one of the typical examples where the cluster and shell-model-like aspects coexist. Since 12 C is the double 
sub-shell closure of the p 3 / 2 shell in the j-j coupling picture, the ground state should contain such the shell- model-like 
I | nature. On the other hand, 3a-cluster structure is also favored in 12 C and is considered to develop in excited states. 

■ Shell model calculations succeeded to reproduce experimental data of many levels in light nuclear region(for example, 
' [1>2]). However, a number of states in light nuclei have been left unsolved in those shell model studies because it is 

difficult for the shell models to describe well-developed cluster states in general. In fact, even the large basis shell 
model calculations failed to describe such the excited states as the 0^" state at E x = 7.65 MeV and the 3j~ state at 
E x = 9.64 MeV near the 3a threshold energy [3]. These states have been considered to be well-developed cluster 
states, and their structures have been studied with 3a cluster models for a long time [4-8]. Recently, Tohsaki et al. 
proposed a new interpretation of the 0^" state as a dilute cluster gas state of weakly interacting 3 a particles [9] , and 

■ discussed it in the relation with the Bose Einstein Condensed(BEC) phase in dilute nuclear matter [10]. Moreover, 
\ there are many experimental reports on the resonances in the excitaion energy region around E x = 10 MeV recently 

[11-14]. Especially, the experimentally suggested + and 2 + states in the broad resonances at E x ~ 10 MeV attract 
an interest because these states are the candidates of the well-developed cluster states. Although the structure of 
these resonances [15-18] has been theoretically studied based on the 3a cluster model approaches, their properties 
have not been clarified yet. The cluster model approaches are useful to describe the details of inter-cluster motion in 
J> \ the 3a system, however, they failed to reproduce such the data as the excitation energy of the 2f and the strengths of 
(3 decays, because these properties are sensitive to such the shell-model-like features as the sub-shell closure and the 
a breaking of 12 C but the 3a-cluster models is not suitable for describing the j-j coupling features nor dissociation 
of the a clusters. 

These facts indicate that the cluster and shell-model-like features certainly coexist in the 12 C system. It is expected 
that the shell-model-like aspects in the low-lying states may have some influence on the excited cluster states through 
the orthogonality and mixing with the lower states. Therefore, it is important to take into account both aspects in the 
systematic study of 12 C. However, there are a few theoretical studies of 12 C where both the cluster and the shell-model- 
like features are taken into account. The method of antisymmetrized molecular dynamics(AMD) has been proved to be 
a powerful approach to describe the cluster and shell-model-like features in general nuclei. Since all the single-nucleon 
wave functions are treated as independent Gaussian wave packets in this method, the existence of clusters are not 
assumed in the model. Instead, if the system favors a cluster structure, the cluster state automatically appears in the 
energy variation in this framework. For the excited states, the method of variation after spin-parity projection (VAP) 
in the framework of AMD is a useful approach [19-21]. The author has performed the first VAP calculations of AMD 
and briefly reported its application to 12 C [19]. Itagaki et al. also have investigated the cluster and shell competition 
in the 12 C system while incorporating the dissociation of a clusters [22]. However, the work was limited to the 
low-lying states. Neff et al. have done the VAP calculations within the fermionic molecular dynamics(FMD) [23], 
which is similar approach to the AMD. Although the FMD wave function is suitable for describing both the cluster 
and shell- model- like structures in principle as well as the AMD, the existence of 3 a-clusters were apriori assumed in 
the practical application to 12 C in [23]. 
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In the present work, we studied the structure of the ground and excited states of 12 C with the method of the VAP in 
the AMD framework. The method is almost the same as that applied in our previous work [19]. We calculated various 
quantities such as radii, transition strengths and densities, and compared them with the experimental data and other 
theoretical calculations. The structures of the ground and excited states were analyzed while focusing on the cluster 
structure and the dissociation of the a clusters. We also estimated the partial decay widths of the resonances with 
the method of reduced with amplitudes. 

This paper is organized as follows. In II, we give a brief explanation of the formulation of the present work. 
The adopted effective nuclear interactions are described in III. In IV, we show the calculated results concerning 
such observables as the energy levels, radii and (3 decays as well as the EO and E2 transitions compared with the 
experimental data. In an analysis of the obtained wave functions in V), the shell- model-like and cluster features are 
discussed. In the last section(VI), we give a summary. 



II. FORMULATION 



The detailed formulation of the AMD method for nuclear structure studies is described in [19,20,24-26]. In partic- 
ular, the formulation of the present calculations is basically the same as that described in [19-21]. 
An AMD wave function is a Slater determinant of Gaussian wave packets; 

$amd(Z) = -^=A{ipi,if2, <pa}, (1) 

where A is the antisymmetrizer and A is the mass number. The i-th single-particle wave function is written by a 
product of spatial(0), intrinsic spin(x) and isospin(r) wave functions as, 

= </>x,XiTi, (2) 

4> Xi (rj) cx exp{-!/(rj - ^=) 2 }, (3) 
\ v 

Xi = {\+Zi)Xl + {\-Zi)Xl- (4) 

The spatial part is represented by complex variational parameters, X H , X 2 j, X 3i , which indicate the center of the 
Gaussian wave packet. The orientation of the intrinsic spin is expressed by a complex variational parameter £i, and 
the isospin function is fixed to be up(proton) or down (neutron). We used a common value of the width parameter u, 
which is chosen to be the optimum value of each nucleus. Accordingly, an AMD wave function is expressed by a set 
of variational parameters, Z = {Xi, X2, • • • , X^, £1, £2, • • • , £,a}- 

We performed the energy variation after spin-parity projection (VAP) in the AMD model space. In order to ob- 
tain the wave function for the lowest J 71 " state, we varied the parameters Xj and = 1 <~ ^4) to minimize the 
energy expectation value of the Hamiltonian, ($| Jf |$)/($|$), for the spin-parity projected AMD wave function; 
$ = P//^-,$amd(Z). P'mk' i s the spin-parity projection operator. By the VAP calculation for the J 71 ' eigen state 
with an appropriate K' , we first obtained a set of parameters Z = Zf 71 " for the lowest J n state which is expressed 
by -f > Mi<-'^AMD(Z'/ 7r ). In order to search for the parameters Z for the higher state, the VAP is performed in the 
orthgonal space to the lower states. The parameters Z^ for the nth state are provided by varying Z so as to 
minimize the energy of the wave function; 

m = |iW'®AMD(Z)) - 2^ I^W$AMD(Z fc )) . — — g, (7 j v sy ( 5 ) 

k=1 \ MK' AMD )\ r MK' 9 AMD{^ k )) 

which is orthogonalized to the lower states. 

After the VAP calculation of the states for various J, n and ir — ±, we obtained the optimum intrinsic states, 
^ > AMD(Z^ 7r ), which approximately describe the corresponding J™ states. In order to improve the wave functions, we 
superposed all the obtained AMD wave functions. Namely, we determined the final wave functions for the J£ states 
by simultaneously diagonalizing the Hamiltonian matrix, (P^ Tr R -,<I)AMD(Z^ , . 7ri )|iJ|P^ : ,, ( f>AMD(Z^ 7r3 )), and the norm 

matrix, (P^ K ,$A M v(Z , ^ i )\Pj^ T K ,,<fr AM v(Z'n j 7T: ')), with respect to for all of the obtained intrinsic states and 

(K',K"). Consequently, the J% state are written as, 

\JZ) - E < J n > i. K)\P^ AM u(Z J k ^)), (6) 

i,K 
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where the coefficients c{J^^i^K) arc determined by the diagonalization of the Hamiltonian and norm matrices. In 
comparison of theoretical values with the experimental data, we calculated the expectation values for the corresponding 
operators by the | J*)- 

III. INTERACTIONS AND PARAMETERS 

The effective nuclear interaction adopted in the present work consists of the central force, the spin-orbit force and 
the Coulomb force. We adopted MV1 force [27] as the central force, which contains a zero-range three-body force in 
addition to the two-body interaction. Concerning the spin-orbit force, the same form of the two-range Gaussian as 
the G3RS force [28] is adopted. The adopted interaction parameters are the same as those used in the previous work 
[19]. Namely, the Bartlett, Heisenberg and Majorana parameters in the MV1 force are chosen to be b = h = and 
m = 0.62, and the strengths of the spin-orbit force are taken to be ui = —Uu = 3000 MeV. 

We used the width parameter v = 0.19 fm -2 for the single-particle Gaussian wave packets in Eq. 3 so as to minimize 
the energy of the lowest + state. 

In the previous work [19], we situated an artificial barrier potential in the procedure of the energy variation to 
confine nucleons within the inner region in the same way as Ref. [20]. Since we are interested in the dilute cluster 
states, we performed the energy variation without the barrier potential. We also did the variation with the barrier 
potential and combined the obtained wave functions as the base wave functions in the diagonalization of Hamiltonian 
and norm matrices. 



IV. RESULTS 

The wave functions of the ground and excited states of 12 C were calculated based on the VAP calculations in the 
framework of AMD. In this section, we show theoretical results and compare them with experimental data. 

The wave functions for the lowest ,V states were obtained by the VAP calculation of Pm K ,$amd by using ( J 71 , K')= 
(0+,0), (2+,0), (4+,0), (6+,0), (1",1), (2-,l), (3-, 3), (4", 3), (5", 3), (l+,0). After obtaining the lowest states(Jf), 
we calculated the second and third J 77 states (0^, 0^, 2^, 2^, 4^) with the VAP calculation in the model space 
orthogonal to the obtained lower J 71 states. 

Those VAP calculations were performed without the artificial barrier potential. As shown later, thus obtained 
intrinsic wave functions for the 0^~, 2 J, 2^~, 4^, 6^, l^f, 2^ states are not compact states but spatially expanded states 
where a particles distribute far away from the center. For these states, we also performed the VAP calculations with 
the artificial barrier potential [20] in the same way as the previous work [19] to keep the particles in an intermediate 
distance region. 

The obtained AMD wave functions(<i> AMD (Z^ 7r )) are considered to approximately describe the intrinsic states of the 
corresponding J™ states. The final wave functions of the J™ states were determined by superposing the spin-parity 
eigen states projected from these obtained AMD wave functions so as to simultaneously diagonalize the Hamiltonian 
and the norm matrices. For the 1 + state, we found that the obtained wave function, P^ f l 1 <I , AMD (Z 1+ ), contains a 
significant isospin breaking component. Therefore, we added the mirror wave function Pp^ri^Mi < E > AMD(Z 1+ ) of the 
original one to the set of base AMD wave functions in the diagonalization. Here P p ^ n is the proton-neutron exchange 
operator. As a result, the number i of the superposed AMD wave functions P^^^amd^;^) is 23 in the present 
calculations. 



A. Energies 

The theoretical binding energy of 12 C is 88.0 MeV. Although the calculation slightly underestimates the experimental 
value 92.16 MeV, it can be improved by changing the Majorana parameter m of the interaction. 

In Fig. 1, the energy levels of the present results(AMD) are shown with the experimental data and the other 
calculations of 3a cluster models; the resonating group method(RGM) [4] with the Volkov No. 2 force and the generator 
coordinate method(GCM) [6] with the Volkov No.l force. Recently, further theoretical works were performed based 
on 3a cluster models, where resonant 3a states were treated in more details [15-18]. The energy spectra obtained in 
these works are almost consistent with either of the 3aRGM [4] and 3aGCM [6] results shown in Fig. 1, except for 
the new broad 0+ state predicted in [16]. 

In the ground band, the AMD calculation reproduces well the level spacing between 0+ and 2+ , while all the 3a 
calculations underestimate it. The success of the AMD is owing to the energy gain of the spin-orbit force in the 0+ 
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FIG. 1. Energy levels of the ground and excited states of 12 C. The theoretical levels of the 3qRGM and 3aGCM are the 
results of [4] and [6], respectively. 

with the breaking of the 3a clustering. In other words, this large level spacing is one of the characteristics of the p 3 / 2 
sub-shell closure, which can not be described within the 3a cluster models. 

For the positive parity states with natural spins above the 3a threshold energy (7. 272 MeV), the 0^, , 2%, 2g~, 
4+, 4 2 and 6± states are obtained in the AMD results. The level structure of these states is similar to that of the 
3aGCM calculations [6]. Since the highly excited states, Og , 2^ ,2^ and 6^, have the dominant component with 
an a-particle far from the other two a's, the stability of these states should be analyzed by treating the boundary 
condition of the resonant states carefully. 

The level spectra of the positive parity states around E x = 10 <~ 11 MeV and their properties have not been clarified 
experimentally and theoretically. In the AMD results, the 0^ state appears at 3 MeV higher energy than the 0^" state. 
The J state was predicted also in the 3aGCM calculation [6] , while Kurokawa et al. proposed a new broad + state 
below the + in the recent work with the method of analytic continuation in the coupling constant combined with 
the complex scaling method(ACCC+CSM) [16]. As discussed later, the O3 in the present result corresponds to the 
O3 state of the 3aGCM calculation. The detailed assignment of the excited states is discussed in the next section. 

We obtain the if and \\ states. These states correspond to the 1+(T = 0) at 12.7 MeV and the 1+(T = 1) 
at 15.1 MeV. Since the spin-parity 1+ is unnatural in a 3a system, it is difficult to describe these states within 3a 
cluster models. In fact, the present 1 + states have the non-3a component. The excitation energies of the 1+ states 
are overestimated in the present results. They can be improved by tuning the strength of Bartlett and Hciscnbcrg 
terms. For example, by using the interaction parameters, b = h = 0.2 and m = 0.62, we obtained the excitation 
energy E X {1+) = 13.9 MeV. 

In the negative parity states, we obtained {37/, 4]", 57} and {17/, 27/, 3.7, 47, 57/}, which construct the rotational 
bands, K = 3~ and K = 1~, respectively. These bands are consistent with the 3aGMC calculations. 

B. Radii 

The theoretical values of the root mean square radius(i? m ) of mass distributions are shown in Table I. The calculated 
radius of the ground state is R m = 2.53 fm, which is slightly larger than other 3a calculations and is large compared 
with the observed value R m = 2.32 — 2.33 fm deduced from the charge radius. The radii(i? m ) of the excited states in 
the present results seem to be qualitatively consistent with the those of the 3aGCM calculations. 

Recently, the structure of the 0^" state attracts a great interest, because Tohsaki et al. [9] proposed a new interpre- 
tation of the 0^~ as a dilute a-cluster gas, where weakly interacting 3 a's form an a condensate state. In their work, 
the O2 has the extremely large radius because of its dilute 3a gas-like structure. Although the radius of the 0^ is 
somewhat smaller in the present result than the theoretical values of the 3a calculations, it is still remarkably large 
compared with that of the ground state. In the present framework, the long tail of the inter-cluster wave functions at a 
dilute density may be underestimated because the base AMD wave functions is limited to only 23 Slater determinants 
obtained by the VAP. We think that the dissociation of the a clusters can be another reason for the smaller radius 



TABLE I. Root mean square radius for mass distributions. The theoretical values of the 3a-cluster calculations [4,6,17] are 
also listed. The values for the 3a condensate wave functions in Ref. [17] are those obtained with the Volkov No. 2 force. The 
mass radii(i? m ) for the 3qRGM calculations [4] are deduced from the charge radii(_R c ) with the relation R? c = — Rp, where 
R p is the proton charge radius 0.813 fm. The observed charge radius of the ground state is R c — 2.46 — 2.47 fm [29], which 
corresponds to R m = 2.32 — 2.33 fm. 
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of the Of than the 3a calculations. One of the origins of the dilute character of the Of state is the orthogonality to 
the ground state. In the 3a cluster model, the ground state has the compact 3a wave function, the wave function of 
the Of state tends to avoid overlapping with the compact 3a state to keep the othogonality. However, in the present 
results, the ground state contains the a breaking component as naturally expected in the p 3 / 2 sub-shell closed nucleus. 
Due to the non-3a component, the compact 3a wave function is partially allowed to mix into the Of state, and it 
may reduce the size of the Of state. 

In Fig. 2, the matter density distributions are shown. The Of , 2f and 4f in the ground band have compact density 
distributions. The density distributions of the 1 + states are also compact and its radius is small because of the 
shell-model-like structure with the dominant Ohco components as well as the ground state. In the Of state, the density 
is suppressed in the small r region and spreads to the outer region. The shapes of the densities in the Of , 2 f and 
4f are similar to each other. They have long tails due to the 3a linear-like structure. In the 3]~ and 4j~ states, the 
density at the center r = fm is small, and the positions of the maximum density are large as r <~ 1.5 fm compared 
with those in other states because of a developed triangle 3a cluster structure. 



C. E2 and E0 transitions 



The E2 and E0 transition strengths are listed in Table II. The present B(E2:2^ — > Of) value for the transitions 
in the ground band well agrees with the experimental data and is consistent with other theoretical calculations. For 
the transition between the 2f and the Of , the 3a cluster calculations [6,4] usually underestimate the observed data 
by factor 2^4. In the present results, the B(E2; 2f — > Of ) value is larger than those of the 3a calculations. It 
is because the spatial extension of the 3a distribution in the Of is not so remarkable as the 3a calculations due 
to the a breaking component of the ground state, and can have a significant transition overlap with the 2f state. 
It means that the a breaking component affect the B(E2;2f — ► Of ) as argued in the work with the a + Ap + An 
model by Itagaki et al. [22]. Considering that the present result overestimates the observed B(E2;2f — > Of), the a 
breaking effect might be somewhat too strong. In the E2 transitions from the highly excited states such as the 2f , 
Of, 4f , we obtained remarkably large B(E2) values, B(E2; 2+ -> Of ) = 100 c 2 fm 4 , B(E2; 2+ -> Of ) = 310 e 2 fm 4 and 
B{E2-.A\ — > 2f ) = 600 c 2 fm 4 , due to the large intrinsic deformation of the developed 3a cluster states. Since the E2 
transitions 2f — > Of and 4f — > 2f are especially strong, we consider that those are the intra-band transitions among 
the band members {Of, 2f and 4f }. This assignment of the band members is also consistent with the similarity of 
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TABLE II. Strengths of E2 and EO transitions. The experimental data are taken from Ref. [32]. The theoretical results of 
3aGCM [6], 3qRGM [4] and a + 4p + 4fc calculations [22] are also listed. 
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the density form among these three states(see Fig. 2). The detailed discussion of the band structure is given in the 
next section. 

( X— 2") 

The E2 transition densities Pj~j. (r) are shown in Fig. 3. The density is normalized as, 

B(E2; Ji ^J f ) = ^jL- J rWj=j? (r). (7) 

The shapes of the transition densities for 2f — > Of , 2 J — > Of , and 2f — > Of in the present results are consistent with 
those of the 3aRGM calculations [4]. The transition density for 2f — > Of has a significant strength in the region 
r = 2 ~ 3 fm in the present results, while, in the 3aRGM case, it has a node which suppresses B(E2). 

The elastic form factor for the ground state is shown in Fig. 4 (a). The dip position shifts toward the small q region 
compared with the observed data, because the present calculation overestimates the radius of the ground state as 
shown in Table I. The inelastic form factors for the electric monopole transitions, Of — ► Of and Of — > Of are shown in 
Fig. 4(b). The first peak of the form factor for Of — ► Of well agrees with the observed data for the inelastic scattering 
into the Of state at 7.65 MeV. The second peak at the high q region is smaller than the observed data because the 
calculation underestimates the observed elastic form factor of the ground state in this region. Recently, Funaki et al. 
discussed the dependence of the inelastic form factor on the nuclear size of the Of within the 3a cluster model [30]. 
They argued that the maximum value of the form factor depends on the size. If the Of state with the 3a structure 
has a small size, the inelastic form factor becomes large because of the large transition overlap with the compact 
ground state. In their analysis the observed form factor is well reproduced by the a condensate wave function with 
the radius R m = 3.8 fm. Also in the 3aGCM and 3aRGM calculations, the form factor is well reproduced by the 
Of with R m — 3.5 fm. On the other hand, although the size of the Of is as small as 3.3 fm in the present results, 
the maximum peak height of the calculations is in good agreement with that of the observed inelastic form factor. In 
the present case, the a breaking component in the ground state suppresses the transition into the 3a cluster state, 
because the non-3a component has a small transition overlap with the 3a cluster state. It means that the absolute 
value of the inelastic form factor is not sensitive only to the size of the Of state, but also depends on the a breaking 
component in the ground state. 

The calculated form factor for the Of — > Of has a similar q dependence to that of the Of — ► Of . Its magnitude is 
one order smaller that for the Of — > Of in the present results. It corresponds to the 3 times smaller transition matrix 
M(M);0f -> Of) = 2.0 fm 2 than the M(M);0f -> Of ) = 6.0 fm 2 . The isoscalar monopole strengths to the Of and 
the Of were studied by inelastic 6 Li and a scattering [11,31]. From the isoscalar energy weighted sum rule strength in 
those studies, the ratio B(E0; Of — > Of )/B(E0; Of — * Of) is deduced to be 0.3 ~ 1 by assuming the mirror symmetry. 
The present form factor for the Of — > Of is smaller than these experimental values. In the 3aGCM calculations [6], 
the inelastic form factor for Of — > Of is as large as that for Of — > Of and is consistent with the observed data. In 
order to study the detailed structure of the Of state, more precise experimental data are required. 
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FIG. 3. Calculated E2 transition densities p*j~j}{r) between the + and 2 + states. 
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FIG. 4. Squared (a)elastic and (b)inelastic form factors of the electron scattering from 12 C given by the calculated transition 
densities and by the observation [33] . 
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TABLE III. Calculated log ft values of the Gamov- Teller transitions for f3 + decays compared with the experimental data 
for 12 N(/3+) 12 C* taken from Ref. [32]. 
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D. GT and Ml transitions 

The /3 decay strengths from neighboring nuclei, 12 B and 12 N to the ground and excited states of 12 C are experimen- 
tally known. The observed values are of great help to investigate the structure of 12 C as well as the electric transition 
strengths. However, within the 3a cluster models, it is difficult to discuss the [3 transitions because Gamov- Teller(GT) 
transition matrix for the ideal a cluster of the (0s) 4 configuration is exactly zero due to the Pauli principle. In other 
words, GT transitions to the 3a cluster component are forbidden, but the transitions are caused by the breaking of 
the a clusters. As a result, the (3 decay strengths are sensitive to the a breaking component. 

In Table III, the calculated GT transition strengths for the /3 + -decay from 12 N are shown and compared with the 
experimental data. The ground state(J 7r = 1 + ) of the parent nucleus 12 N is obtained with a variational calculation 
after spin-parity projection to the 1 + with the width parameter v = 0.19 fm~ 2 . The calculated strengths are in good 
agreement with the experimental data. The strong (3 decay to the ground state is because of the rather large mixing 
of the a breaking component, which reflects the shell-model- like feature of the j-j coupling. The GT transitions 
to the Of and O3 states have significant strengths though these states are dominated by the developed 3a cluster 
states. The strengths of the transitions to these excited + states come from the mixing of the wave functions in 
the diagonalization. Before the diagonalization, it was found that the GT transitions are weak for the original single 
AMD wave functions, P^ K <&amd(Z 2 ~3) for the 0^~ and Of , where the 3a cluster structure is developed. After the 
diagonalization, the P^-$amd(ZJ + ), P^k ^amd(Z2 + ) and P^ K ^p L - MD (Z^ + ) mix in the superposed wave functions, 
and the final wave functions \0f) and \0f) contain small fraction of the Pj^-<I>AMD(Zi + ). As a result, the a breaking 
component, which is originally included in the ground state wave function P^-$amd(Z5 + ), is contained in the \0f) 
and the \0f) through the superposition of the wave functions. This is the origin of the significant GT transitions to 
the Of and the Of. The GT transition to the 2f is weaker than that to the Of. It is reasonable because the 3a core 
structure enhances in the 2f compared with the 0^~ as shown in the next section. The strong GT transitions to the 
lf{T = 0) state at 12.71 MeV is described well by the present calculations reflecting the intrinsic spin excitation in 
this state. 

The calculated Ml transition strengths are shown in Table IV compared with the experimental data. The results 
for the transitions from 1+ and if are in good agreements with the experimental data for the lf(T = 0) at 12.1 MeV 
and the lf(T = 1) at 15.1 MeV. Therefore, we assigned the present if and if states to the observed Vf(T = 0) and 
lf(T = 1) states, respectively. 

V. DISCUSSIONS 
A. Intrinsic structure 

The present wave functions are given by the linear combination of the spin-parity eigen wave functions projected from 
23 AMD wave functions {<&amd(Z^ 71 ')}. Here, $amd(Z^ 71 ') is the optimum AMD wave function, which is obtained by 
the VAP calculation for the J% state. Since each AMD wave function is written by a single Slater determinant, it is easy 
to analyze its intrinsic structure. Even though those bases are superposed by the diagonalization of the Hamiltonian 
and norm matrices, the ^AMDfZ^) is mostly the dominant component of the J™ state after the diagonalization. 
Then we firstly discuss the intrinsic structure of the <j>AMD(Z^). In Figure 5, we show the density distributions of the 
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TABLE IV. Ml transition strengths. The experimental data are taken from Ref. [32]. The unit is the nuclear unit, 
(eh/2Mc) 2 . The experimental excitation energies of the initial and final states are shown in parenthesis. 



B(M1) 

J? J] Exp. Theory(AMD) 

if (12.1 Me V) 0f(0 MeV) 0.015 (0.002) 0.011 

if (12.1 MeV) 2+(4.4 MeV) 0.0081 (0.0014) 0.008 

L+(15.1 MeV) 0+(0 MeV) 0.95 (0.02) 0.17 

1+(15.1 MeV) 2+(4.4 MeV) 0.068 (0.009) 0.09 

1+(15.1 MeV) 0j(7.7MeV) 0.23 (0.04) 0.014 

1+(15.1 MeV) 1+(12.1 MeV) 3.6 (1.1) 2.5 



intrinsic wave functions, < I>AMD(Z^ 7r ). As seen in Fig. 5 (al) and (e3), the Of and the 1 + states have compact density 
distributions and no developed cluster structure. In the 2f and the 4f states of the ground band built on the Of 
state, the 3a cluster core appears(Figs. 5-bl and 5-cl). These states are considered to be the SU(3) limit 3a cluster 
states because the spatial development of the clustering is not remarkable. In the excited states, a variety of spatial 
configurations of the developed 3a cluster structure appears. The 3a in the 0;f (Fig. 5-a2) has an isosceles triangle 
configuration which is close to an equilateral triangle. The similar configurations of the 3a are seen also in the 3f , 4f 
and 5f states (Fig. 5-d3,5-el,5-c2). In the 0^ state(Fig. 5-a3), the 3a shows a rather linear-like configuration, where 
the largest angle of vertices is larger than 120 degree. We found such the linear- like structure also in the 2 J, 6f , If 
and 27". In these states, a a cluster seems almost escaping. In Fig. 5(fl)-(g3), we show the intrinsic structure obtained 
by the VAP with the artificial barrier potential as explained in II. Although the a particles are more confined than 
the case without barrier, the developed 3a cluster structures are still formed in these wave functions. 

The spatial 3a configurations of the 3~ and the l~ states in the present results are consistent with the energy 
minimum states in the isosceles configurations of the 3a model space [6] . On the other hand, the present configuration 
of the O2 state does not correspond to the energy minimum but seems to be the second minimum state in the 3a 
models space for the 0^" [6]. After the superposition of the basis, the features of the |0^) are similar to those of 
the 3aGCM calculations, because various configurations of the 3a cluster states mix in both the present results and 
3aGCM calculations. 

As explained in Sec. II, we describe the final wave functions | J*} with a linear combination of the obtained wave 
functions, P^ f 1 ^> < f ) AMD(Z / f 171 ), each of which is the spin-parity eigen wave functions projected from a AMD wave 
function. The coefficients of Eq. 6 are determined by the diagonalization of the Hamiltonian and norm matrices 
with respect to all the obtained AMD states. In Table V, the amplitudes of the dominant components in the final 
wave functions \J£) are listed. In many cases, the dominant component in the \J£) is the P^-^amdIZ^ 7 ^) with 
J7£. = J%, and therefore, the $amd(Z^) is regarded as the approximate intrinsic state of the state. 

The \0f), \2f) and |4+) in the ground band contain about 90% component of the dominant P^-$AMD(Z^' r ). In 
the 3a cluster states, the amplitudes of the dominant component decrease mainly due to the mixing of various 3a 
configurations. Especially, in the |0j), the amplitude of the Pqq~<I>amd(Z2 + ) is only 0.49, and it contains other 3a 
configurations such as P 0[ f $amd(Z3 + ). The mixing of the different 3a configurations in the jOj } enhances the loosely 
binding nature of the 3a particles, and increases the ,5- wave component of the 8 Be(0 + )-a motion. On the other hand, 
the dominant amplitude in the |3f ) is 0.82, which is large compared with other 3a cluster states. It may indicate 
that this state still has a nature of the rigid structure of the triangular 3a configuration. The intrinsic state of the 
dominant component in the |4j) and |6f ) is $amd(Z 1 + ), where an a particle is almost escaping and locates at the 
distance r ~ 7 fm. Their stability against a decays should be carefully discussed to confirm the possible existence of 
these states. Since the |2^) is dominated by P^-'&amdIZ^ 7 ^) with JTJj.* = 2 3" and K = 2, this state is considered 
to be the band head state of a side-band K n = 2+. In the negative parity states, the intrinsic structures of the 
*AMD(Zf)'B with J£ — 3f , 4 1 i 5-l are similar to each other, and they construct the — 3 band; {3 1 , 4 X , 5 2 }. 
On the other hand, the $amd(Z^ 7T )'s for J£ = If and J£ = 2f have almost the same linear-like structures and form 
the K n = 1~ band, { If, 2f , 3f , 4f , 5f }. These negative parity bands, K n = 3~ and K 71 = 1~ were suggested also 
in 3aGCM calculations. 

In Table VI, we show the expectation values of the total spin(J), the intrinsic spin(S), the angluar-momcntum(L), 
and the principal quata of harmonic ocsillator(A rho ) for the spin-parity projected states P^-<I>AMD(Z^ 7r ). Here, the 
ocsillator quanta A^ ho is defined as, 
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FIG. 5. Mass density distributions of the intrinsic state $amd(Z^"") before the projection for the states obtained by the 
VAP calculations with respect to the J*. Figures (fl)-(g3) are the mass density distributions of the states obtained with the 
artificial potential barrier. The intrinsic system is projected onto a plane which contains the major axis of the intrinsic state. 
The density is integrated along a transverse axis perpendicular to the plane. The size of the frame box is 10 fmxlO fm. 
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TABLE V. Overlap between Ptf K $ A md(Z^) and the final wave functions \J„). The \JZ) is obtained by linear combination 
°f Pmk $AMD(Zj^ lr< ) as explained in Eq. 6. K is taken to be the same value as chosen in the VAP calculations. The 
overlap with the largest component in the |J£) is also shown in case that the overlap with other components is larger than 
I(WmV1>amd(Z^))| 2 . 
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(8) 

where the width parameter v is chosen to be the same value as that of the single-particle Gaussian wave packets. We 
also show the expectation values of those operators for the final wave functions \J„) after the diagonalization in Table 
VII. 

Non-zero values of the intrinsic spins indicate the breaking of 3a clustering. Due to the 3a cluster structure, the 
(S 2 ) values are almost zero in many states before the diagonalizaton except for the + and 1 + states as seen in 
Table VI. The origin of the non-zero (S 2 ) in the + state is the mixing of the p 3 / 2 sub-shell closure component. In case 
of the 1 + state, non-zero (S 2 ) originates in the intrinsic spin excitation. More detailed discussions of the a-cluster 
breaking is given later. 

The proton total spin (J 2 ) and neutron total spin (J 2 ) are the smallest in the + state. It is because of the significant 
P3/2 sub-shell closure component. On the other hand, they are quite large in the 0^" and the Ojj" states compared with 
the + . This comes from the component of the orbital-angular momenta L p 7^ and L n 7^ 0, which are caused by the 
strong correlations of nucleons in the developed a clusters. In case of the + state, those values becomes larger in the 
|0^) after the diagonalization of the basis than those in the P u0 f $AMD(Z2 + ) before the diagonalization. It indicates 
that the further developed 3a configurations mix in the |0^") through the superposition of the basis. 

The oscillator quanta A^ ho is equivalent to the principal quantum number of the spherical harmonic oscillator with 
the oscillation number u> = 2Kv/m(m is the nucleon mass), and its expectation value indicates how large the mixing 
of higher shell components is in the expansion of spherical harmonic oscillator basis. In the ground band members, 
+ , 2 + , 4 + states, the values of the oscillator quanta for protons {(N^ )) and neutrons((iV^ LO )) nearly equal to the 
minimum value 4, which correspond to the OHuj states. On the other hand, those for the spatially developed 3a 
cluster states such as the 0^, O3", 17 and 37 are much larger than 4. These cluster states should be written by linear 
combinations of a large number of the higher shell configurations in the expansion of the harmonic oscillator basis. 
Therefore, it is natural that these cluster states can not be described even by the large-basis shell model calculations 
[3]. 
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TABLE VI. Expectation values of the squared total spin(J 2 ), the squared intrinsic spin(S 2 ), the squared 
angular-momentum(L 2 ), and the principal quanta of harmonic oscillator(iV ho ) for the spin-parity projected states 
Pmk 5amd(Z^ x ) before the diagonalization. K is chosen to be the same value as used in the VAP calculations. The cal- 
culated values for protons(J p , S p , L p , N p °), neutrons(J„, S„, L n , N n °), and the total system(S, L) are shown. 
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TABLE VII. Expectation values of the squared total spin(J 2 ), the squared intrinsic spin(S 2 ), the squared 
angular-momentum(L 2 ), and the principal quanta of harmonic oscillator(iV ho ) for the final wave functions \ JZ) obtained by the 
superposition of the wave functions after the diagonalization. The calculated values for protons(J p , S p , L p , N p °), neutrons(J n , 
S„, L„, Nn°), and the total system(S, L) are shown. 
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B. Breaking of a clusters 



In the present framework, the existence of a clusters is not assumed, but the formation and breaking of the clusters 
are automatically incorporated in the energy variation if such the structures are favored. We can estimate the a 
breaking components with the non-zero values of the squared intrinsic spins, (S 2 ), listed in Table VI and VII. 

Let us first discuss the cluster and non-cluster components in the wave functions by analyzing the results shown 
in Table VI before the diagonalization. The states except for the Of, 2f, 4^ and if have almost zero values of (S 2 ) 
because of the well-developed 3a cluster structures. It means that the a clusters are not broken in these cluster 
states. On the other hand, (S 2 ) = 1.6 in the 0^ state indicates the significant a breaking component. It is natural 
because p 3 / 2 sub-shell closed configuration is favored in the ground state. If the ground state is the pure p 3 / 2 sub-shell 
closed state, the expectation values of the squared total spins for protons((J 2 )) and neutrons ((J 2 )) should be zero. 
Considering the non-zero values((J 2 ) = (J 2 ) = 0.8), the 0^ state is considered to be a mixture of the 5{/(3)-limit 3a 
state and the p 3 / 2 -shell closed state. In other words, the 3a clustering is partially broken due to the spin-orbit force in 
the ground state. The 1 + is the non 3a state with intrinsic spin excitations as seen in its spin magnitude ((S 2 ) = 3.4). 
It is consistent with the fact that the J T = 1 + is unnatural in a 3a system. 

Next we look into the a breaking in the | J%) after the diagonalization. By comparing the values in Table VII with 
those in Table VI, we found that the a breaking component in the original wave function of the ground state mixes 
into the |0^") and IO3") because the (S 2 ) values of the |0^") and | O3") become large due to the diagonalization. In the 
O2 state, the angular momentum (L 2 ^) for protons(neutrons) increases as well as the intrinsic spins (S 2 ) due to the 
superposition of the wave functions. It indicates that not only the a-breaking component but also the component of 
further developed 3a clustering increases in the 0^ state through the diagonalization. 

As mentioned before, the GT transition strengths are good probes to investigate the mixing of cluster and non- 
cluster components, because the transitions to the ideal 3a cluster states are forbidden, and therefore, the strengths 
directly reflect the a breaking component. In the experimental observations, (3 decays from 12 N to the ground state of 
12 C are rather strong as shown in Table III. Even for the transitions to the excited states, Of and O3", the significant 
strengths of the (3 decays were observed. These facts indicate that the a breaking component is contained in the Of 
and Ojj" states as well as the 0^ state. In the present study, we found that the a breaking is significantly included 
in the original 0^ wave function due to the spin-orbit force, while the original 0^ and 0^ wave functions before the 
diagonalization have the well-developed 3a structure but no a breaking. When we calculated the GT transitions for 
the original wave functions before the superposition of the basis, we obtained almost forbidden GT transitions into 
the Of and 0jj~ states because of the 3a cluster structures, and failed to reproduce the experimental log ft values. The 
a breaking components are slightly contained in the 10^ ) and IO3") after the diagonalization through the mixing of 
the original 0^ wave function, which includes the p 3 / 2 sub-shell closed configuration. As a result of the mixing of the 
a breaking components, the log ft values are well reproduced by the present wave functions after the diagonalization. 

The a breaking gives important effects on various properties of the ground and excited states of 12 C as well as 
the GT transition strengths. As mentioned before, the level spacing between the 0+ and 2+ states is as large as the 
experimental data because the ground state gains the spin-orbit force due to the a breaking. This large level spacing 
can not be described by the 3a cluster models. The a breaking component of the ground state should affect also the 
properties of the excited + states. Within the 3a model space, the ground state is almost the (3)-limit 3a state, 
while the excited 2 state is an optimum solution in the model space orthogonal to the 3a ground state. It means that 
the compact 3a state is forbidden for the Of state to satisfy the orthogonality to the ground state. In the present work, 
since the ground state contains the a breaking component, the S?7(3)-limit 3a state is partially allowed for the Of 
state. We consider that this is one of the reasons for the smaller radius of the 0^ state in the present calculations than 
the values obtained by the 3a calculations. The reduction of the Of size due to the non-3a component in the ground 
state should be important also in the description of the experimental value of the E2 strength. The experimental 
data of B{E2; Of — ► 2f) is usually underestimated by the 3a calculations. In the present study, we can describe the 
increase of i?(i?2;0 2 f — > 2+) by the effect of the size reduction, which enlarges the overlap of the 0^ with the 2f. 
Since the present i?(_E2;0 2 f — > 2+) value is larger than the experimental one, the effect of the a breaking might be 
overestimated. 

Let us discuss the a breaking effect on the inelastic form factor to the state. Funaki et al. showed that the 
maximum value of the inelastic form factor is sensitive to the spatial extension of the 0^ state. In their analysis 
based on the 3a condensate wave function, the magnitude of the form factor decreases with the enhancement of the 
O2 size. They succeeded to reproduce the experimental data of the inelastic form factor by the 3a condensate wave 
function with a large radius R m = 3.8 fm. On the other hand, as shown in Fig. 4, the present results well reproduce 
the maximum value of the form factor, though the root-mean-square radius of the 0^ is R m = 3.3 fm which is smaller 
than that of their calculations. The E0 strengths should be sensitive also to the structure of the ground state. If 
the ground state contains the a breaking component, the inelastic transition strength to the excited 3a state should 
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TABLE VIII. The calculated reducued a-decay widths(# 2 ), the amplitudes (ryh{r)), S factors, cluster probabilities, a-decay 
widths for the 8 Be(0 + )+a channel. The assumed excitation energies E x for the O3 and the 2f states are shown in parenthesis. 
The calculated a decay widths are obtained with the experimental energy E a of the relative motion between 8 Be(0 + ) and a. 
The experimental widths [32] are also shown. 
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decrease, because the non-3a component reduces the transition overlap between the ground state and the a cluster 
states in general. This is one of the reasons for the suppression of the inelastic transitions in the present calculations, 
which reproduce the magnitude of the inelastic form factor for the — > Of. 

It should be stressed that the a breaking component in the ground state is directly reflected in the inelastic transition 
strengths from the ground state, and also can have an influence on the structure of the excited + states. 

C. 3a cluster features 

In many theoretical works, the 3a cluster features of 12 C have been discussed by many groups for a long time. 
Especially, the + and the 2 + states above the 3a threshold energy have attracted great interests. Recently, there 
are some experimental reports on the + and 2 + around E x = 10 MeV [11—14]. However, the level structure and the 
assignment of these resonances have not been clarified yet. In the present work, various 3a cluster structures appear 
in the excited states of 12 C. We, here, analyze the structure of the + , the 2 + and the 4 + states by extracting the 
8 Be(0 + )+a components and estimate the a decay widths with the method of reduced width amplitudes. We give a 
discussion of the level assignment and the band structure later in VD. 

We extracted the reduced width amplitudes for the 8 Be(0 + )+a, and calculated the spectroscopic factors(S) and 
cluster probabilities. For simplicity, we assumed the (0s) 4 wave function for the a particle and (3)-limit cluster 
wave function for the 8 Bc(0 + ) by using the same width parameter v as the present AMD wave function of the 12 C. 
The definitions and the detailed method of the calculations are explained in Rcf. [21]. 

Figure 6 shows the reduced width amplitudes ryz,(r) of the + , 2+ and 4+ states. The calculated S factors and 
cluster probabilities are listed in Table VIII. In the Of, 2f and 4f states, the node numbers of the cluster wave 
functions are 2, 1 and 0, respectively, which correspond to the lowest allowed node numbers for the 8 Be-a relative 
motion. In the higher J n states, the node number increases one by one and the position of the largest peak shifts 
toward the large distance region. The S factors in Table VIII indicate the spatial developments of the 8 Be(0 + )+a 
clustering. The S factor is the largest in the Of state. Among three 2 + states, S is the largest in the 2f state. In the 
Of state, the S factor is not so large even though this state is the developed 3a cluster state. It is because the main 
component of the 0^" is a rather linear-like 3a, which results in the reduction of the 8 Bc(0 + )+a component because 
of the large mixing of the 8 Be(2+)+a with D -wave relative motion. These behavior of the cluster wave functions 
rj/z^r) for the Of, Of, Of, 2f, 2f and 2f corresponds well to the 3aGCM calculations by Uegaki et al. [6], though 
the present amplitudes are slightly smaller than those of the 3aGCM results. Therefore, the basic features of the 3a 
cluster components in the present results are considered to be similar to those of the 3aGCM results. 

The calculated reduced widths(# 2 (a)) are shown in Table VIII. The large values of 0„(a) in the Of, Of and 2f 
states imply the spatial development of the a cluster in these states. We estimated the 8 Be(0 + )+a decay width with 
the method of reduced width amplitudes in the same way as in Ref. [21]. In the present estimation, we temporary 
assigned the Of and the 2f to the observed states around E x = 10 <~ HMeV and chose the energy of the relative 
motion E a to be 3 MeV. The calculated a decay widths for the O3 and the 2f are 400 keV and 460 keV, respectively. 
In these energy region, the partial width for the 8 Bc(0 + )+a decay is expected to dominate the total a decay widths 
as shown in Ref. [6]. 
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FIG. 6. a-reduced width amplitudes ry(r) of the + , 2 + and 4 + states for the [ 8 Be(0 + ) ® = J) channel. 



17 



D. Assignments of + and 2 + resonances 

The + and 2+ states around E x = 10 McV have been experimentally and theoretically studied for a long time. 
We here consider the correspondence of the present results with other theoretical calculations and the experimental 
reports. 

There exists the broad resonance at E x =10.3 MeV with the width r ~ 3 MeV, which has been ambiguously 
assigned to a + state [32]. Because of the existence of this broad resonance, there remain uncertainties on the 
spectroscopy in this energy region. In the inelastic 6 Li and a scattering [11,31], a broad + state has been observed 
around E x — 10 MeV. Its energy and width are consistent with the broad + at E x =10.3 McV. The + state around 
10 <~ 11 MeV has been also observed in the recent experiments of the /^-delayed 3a particle decay from 12 B and 12 N 
[12,14]. In the argument in [14], interference between this state and the 0^(7.65 MeV) has so significant influence on 
the spectrum and it makes difficult to determine the resonance energy. Recently, the 2\ state has been discovered at 
E x ~ 10 MeV by the inelastic a scattering [13]. The spectra of this state are overlapped by the broad peak of the + 
at 10.3 MeV, and the fitting parameter for the 2\ is E x = 9.9 ± 0.3 MeV and Y = 1.0 ± 0.3 MeV. 

The + state above the 3a threshold has been long discussed also in many theoretical studies. Morinaga et al. [34] 
proposed a linear-chain 3a structure(LCS) of the 0^" at 7.65 MeV. The possibility of the LCS has been discussed based 
on the analysis of the a decay widths by Suzuki et al. [35] , and it was argued that the LCS is inappropriate to explain 
the large experimental widths of the + (7.65 MeV) and the + (10.3 MeV). In the microscopic and semi-microscopic 
3a cluster model calculations [4,6,15,36], it was found that the 0^" has the developed 3a structure with the dominant 
[ 8 Bc(0 + ) <g> I — 0] j = o and has the nature of the weakly interacting 3a particles. This is quantitatively consistent with 
the features of 3a clustering in the 0^" state of the present results. Recently, the new interpretation of the 0^" as a 
3a-cluster gas was proposed [9,36]. Also in the present study, we found the dilute 3a gas aspects of the 0^", though 
its radius is somewhat smaller than that of the 3a calculations. In these theoretical studies, it is considered that the 
3a structure of the 0^ has been qualitatively established. 

The 0^ state 4 MeV above the 0+ state was predicted in the 3aGCM [6], 3aRGM+CSM [37], and 3aOCM+CSM 
[15] calculations. In the comparison with the experimental data of the + (10.3 MeV,r oxp =3 MeV), the theoretical 
energy is rather higher than the observed data, and the calculated width r ca i ~ 1 McV is narrower than the r eX p — 3 
MeV. Also in the present results, we obtained the 0g state a few MeV above the 0^" state. The structure of this state 
is similar to the 0$ state of the 3aGCM calculations. As mentioned before, the dominant 3a configuration of the 
0^ is rather linear-like as well as the 3aGCM calculations, where the 0^~ state is dominated by [ 8 Be(2 + ) ® I = 2] j = o 
and has less [ 8 Be(0 + ) <g> I = 0]j = o component. This results in the suppression of the a decay width. Comparing with 
the observed broad width, it seems that the theoretical 0% state should not be directly assigned to the experimental 
+ (10.3 MeV,r oxp =3 MeV). Recently, Kurokawa et al. predicted a new broad + resonance below the predicted O3" 
state. It may suggest that another broad resonance might coexist with the predicted narrower O3" state in this energy 
region. The possible coexistence of two + resonances in the broad + spectra at 10.3 MeV with r exp = 3 MeV 
has not been experimentally excluded. By speculating the coexistence of the two + states in this energy region, we 
conjecture that the Ojj" in the present result may correspond to the narrower state of the two + resonances. There 
exists some experimental information concerning the total strengths of the transitions to the broad + (10.3 MeV). 
The inelastic 6 Li and a scattering data imply that the isoscalar monopole transition strengths to the + (10.3 MeV) 
is 30-100% of that to the + (7.65 MeV), while in the present results it is 10% of the strength to the 0^ state. On 
the other hand, the calculated B{GT) for the GT transition from 12 N to the 0^" is half of the experimental strengths 
for the transition to the + (10.3 MeV) (see Table III). In order to clarify the properties of the 0^ state, more precise 
experimental data and further theoretical studies are required. 

Next, we discuss the 2\ state. The cluster models have predicted the 2\ state around E x = 10 MeV as a member 
of the rotational band built on the Oj" state [4,6,15,17]. The recent experimental report on the 2^ at 9.9 ± 0.3 MeV 
with T = 1.0 ± 0.3 [13] is consistent with the theoretical energy and width of the predicted 2 J in Refs. [15,17] and 
also with the present results. The theoretical 2^~ state is dominated by [ 8 Be(0 + ) ® I = 2] j =2 as shown in Table VIII 
as well as in the 3a models [6,15,17]. Therefore, in the weak coupling picture of a 8 Be+a system, the 2^~ state is 
regarded as the rotational member upon the 0^" state, which has the dominant [ 8 Be(0 + ) <g> I — 0]j = q. The E2 strength 
from the 2\ to the 0^" is large as B(E2; 2^ — ► 0^ ) = 100 e 2 fm 4 in the present calculations. However, by analyzing the 
E2 transitions and the intrinsic structure of the + and the 2 + states, we can extract an alternative interpretation 
for the band assignment of the 2^ . Namely, the 2 J and the J can be the rotational band members of the linear- like 
3a structure. In fact, we obtained the remarkably large B(E2) values, B(E2; 2\ — > O3") = 310 c 2 fm 4 . Moreover, the 
dominant intrinsic wave functions of the 2^ and the 0^j~ are quite similar to each other. They show the linear-like 3a 
structure, where the largest angle of vertices of the triangle 3a configuration is larger than 120 degree (Figs. 2-a3 and 
2-b2). The overlaps of the |2^~) and |0^") with the wave function P^ a: $amd(Z§ + ) projected from the single AMD 
wave function (Fig. 2-a3) are more than 60%. As a result of the similar intrinsic structure with the large deformation, 
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B(E2; 2^ — > O3") is extremely large, and the shape of the density p(r) of these two states are quite similar to each 
other(Fig. 2). Consequently, in a collective picture, we consider that the and the 2 J states are regarded as the 
members of the rotational band, which is formed by the deformed intrinsic state with the linear-like 3a structure. We 
expect that the real 2 \ state may have an intermediate rotational feature of this strong coupling aspect built on the 
O3 state and the weak coupling one built on the {)\ state. Unfortunately, the experimental information for the 2^~ 
and the O3 state is too poor to justify the band assignment of the 2% state. 

VI. SUMMARY 

We investigated structure of the ground and excited states of 12 C. The present theoretical method is based on 
the VAP calculations in the framework of AMD, which can describe both cluster and shcll-modcl-likc aspects. The 
present results systematically reproduce the various observed data such as the energy levels, electric and j3 transition 
strengths. It was found that the ground state is the mixture of the shell-model-like state and the 3a cluster state, 
while the developed 3a cluster structures appear in the excited states. We also obtained the 1+ states with non-3a 
cluster structure. 

We discussed the a breaking components. It is important that the a breaking significantly affects not only the 
low-lying states but also the excited 3a cluster states. Firstly, the experimental large level spacing between the 
and 2^ states is described by the energy gain of the spin-orbit force due to the a breaking in the which contains 
the p3/2-shell closure component. It was found that the a breaking component is slightly contained even in the 
excited states with the developed 3a cluster structure through the mixing of the ground state wave functions with 
the p 3 / 2 -shell closure component. As a result of the mixing of the a breaking component, the significant (3 decay 
strengths are well reproduced by the present wave functions. Moreover, it is important that the non-3a component 
in the ground state may change the 3a cluster structure in the excited + states through the orthogonality to the 0+ 
state. 

In the analysis of the intrinsic structure and the 3a clustering in the present results, it was found that the features 
of the 3a cluster components well correspond to those of the 3aGCM calculations [6]. As well as in the 3aGCM, 
we obtained the Oj, 2 J and 0^ states. The 0^, which corresponds to the observed + (7.65 McV), is dominated by 
[ 8 Bc(0 + ) ® I = 0] j = o, and shows the trend of cluster-gas features. The predicted 0^ can not be directly assigned to 
the experimental broad + resonance (10. 3 MeV). We expect that the corresponding + spectra might be overlapped 
by the broad + (10.3 MeV). The predicted 2\ state should be assigned to the recently observed 2 + state at E x = 9.9 
MeV [13]. In the weak coupling picture of 8 Be+a, This 2^ can be regarded as the rotational member built on the 
O2 state, because this state is dominated by [ 8 Bc(0 + ) ® Z = 2] ,7=2. We also proposed an alternative interpretation for 
the band assignment of the the 2^ states based on the remarkable B(E2; 22^^03^) and the similarity of the intrinsic 
structure with the O3 in the present results. Namely, in a collective picture, the O3 and the 2\ states are considered 
to be the members of the rotational band, which is formed by the deformed intrinsic state with the linear-like 3a 
structure. 

In the present work, the number of the base wave functions is limited and the 3a configurations with the large 
distance are not incorporated enough. Therefore, the description of the detailed resonant behavior is insufficient in 
the present framework. Especially, for the broad resonances with 3a clustering, the boundary conditions should be 
carefully treated in the calculation of the excitation energies and the widths. The broad + and 2 + resonances have 
been theoretically studied with the CSM and ACCC methods in the 3a cluster models. Such methods are promising 
in further study on these resonances. 

In conclusion, we should stress that the cluster and the shell-model-like features coexist in 12 C. Since these two kinds 
of nature interplay with each other, one should take into account both the features appropriately for a systematic 
study of 12 C. We expect that it is also significant for clarification of the properties of the 3a resonances. 
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